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ABSTRACT 

We present optical spectroscopic identifications of hard X-ray (5-10 keV) selected sources belonging to 
the HELLAS sample obtained with BeppoSAX down to a 5-10 keV flux limit of /5-iofcey~3 x 10“^'^ erg 
cm“^ s“^. The sample consists of 118 sources. 25 sources have been identified trough correlations with 
catalogues of known sources. 49 have been searched for spectroscopic identification at the telescope. 13 
fields resulted empty down to R=21. 37 sources have been identified with type 1 AGN and 9 with type 
2 AGN. The remaining are: 5 narrow emission line galaxies, 6 Clusters, 2 BL Lac, 1 Radio Galaxy and 
1 Star. Combining these objects with other hard X-ray selected AGNs from ASCA and HEAOl, we 
find that the local luminosity function of type 1 AGN in the 2-10 keV band is fairly well represented 
by a double-power-law-function. There is evidence for significant cosmological evolution according to 
a pure luminosity evolution (PLE) model Lx{z)(x{l -I- z)^, with fc=2.12 and k=2.22 (ct/c~ 0.14) in a 
(Dm,U a)=( 1.0,0.0) and in a (Dm,DA)=(0.3,0.7) cosmology, respectively. The data show an excess of 
faint high redshift type 1 AGN which is well modeled by a luminosity dependent density evolution 
(LDDE), similarly to what observed in the soft X-rays. However, in both cosmologies, the statistic is 
not significant enough to distinguish between the PLE and LDDE models. The fitted models imply a 
contribution of AGNl to the 2-10 keV X-ray background from 35% up to 60%. 

Subject headings: cosmology: observations — galaxies: distances and redshift — galaxies: evolution — 
galaxies: active — quasars: general 


1. INTRODUCTION 

AGN have first been discovered in the radio and soon 
after searched in the optical band. Gonsequently, they 
have been classified using their optical characteristics and 
mainly divided into two categories: type 1 (AGNl) and 2 
(AGN2) according to the presence or not of broad emission 
lines in their optical spectra (we will keep this definition 
of AGNl throughout this paper). 

Before the advent of the last generation of hard X-ray 
telescopes, AGN samples where predominantly based on 
AGNl selected either in the optical or, later on, in the 
soft X-rays by Einstein and ROSAT. In these bands the 
evolution of AGNl has been well measured (see e.g. Della 
Geca et al. 1992; Boyle et al. 2000; Miyaji, Hasinger, & 
Schmidt 2000). On the contrary the production of samples 
of AGN2 has been difficult at any wavelength and limited 
to few local surveys. 

The general picture was in favor of a model in which 
AGNl objects were associated to AGN with low absorp¬ 
tion in the hard X-rays while AGN2 to obscured sources 
with large column densities and spectra strongly depressed 
in the soft X-rays, as expected in the unification models 
(e.g. Antonucci 1993). 

In the last decade the advent of the ASCA and Bep¬ 
poSAX satellites has allowed for the first time the detec¬ 


tion and identification of AGN as the main counterparts of 
hard (2-10 keV) X-ray sources down to fluxes ~ 5 x 10“^^ 
erg cm“^ s“^, more than 2 orders of magnitude fainter 
than HEAOl (Wood et al. 1984). These identifications 
accounted for about 30% of the 2-10 keV hard X-ray back¬ 
ground (Ueda et al. 1998; Fiore et al. 1999). Recently the 
new generation of X-ray satellites such as Chandra and 
XMM-Newton, have reached fluxes 100 times fainter, iden¬ 
tifying hundreds of sources and almost resolving the hard 
(2-10 keV) X-ray background (e.g. Mushotzky et al. 2000; 
Fiore et al. 2000; Giacconi et al. 2001; Hornschemeier et 
al. 2001; Hasinger et al. 2001; Tozzi et al. 2001; Baldi et 
al. 2001). 

Thanks to their excellent angular resolution (^1-5"), 
the first spectroscopic identifications projects have been 
able to observe faint (I~23) optical counterparts. At vari¬ 
ance with the classical type-l/type-2 model in the optical, 
a significant number of the counterparts (~30%) resulted 
to be apparently optical normal galaxies, with X-ray lumi¬ 
nosities — 10"^^ erg s“^ typical of AGN activity, 

and moreover part of the optical type 1 AGNs resulted to 
be absorbed in the hard X-rays (see e.g. Fiore et al. 2000; 
Barger et al. 2001; Tozzi et al. 2001; Hornschemeier et al. 
2001; Gomastri et al. 2002). 

These observations have complicated the picture of the 
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AGN model. In this framework the computation of the 
density of AGN has become an even more difficult task. 
In fact, it is not clear how to classify the sources and to 
take into account the selection biases introduced by the 
observation in the 2-10 keV range, where the absorption 
still play a relevant role. 



Fig. 1.— The flux distribution of the 118 sources from 
HELLAS used in this analysis (continuous line). The 
dashed line is the distribution of the 74 sources for which a 
spectroscopic identification campaign has been carried out 
{empty fields included). The distribution of the 13 sources 
for which we end up with an empty field is represented by 
the hatched histogram. 


These recent deep surveys with Chandra and XMM- 
Newton have reached fluxes ~ 5 x 10“^® erg cm“^ s“^(2-10 
keV) in quite small areas (less than 1 deg^). As a conse¬ 
quence these surveys are not able to provide statistical 
significant samples at brighter fluxes (~ 10“^® erg cm“^ 
s“^; 5-10 keV) where the density of sources is about 5/deg^ 
(Fiore et al. 2001) and tens of square degrees are to be 
covered. Such data are necessary to provide large numbers 
of spectroscopic identified sources in a wide range of X-ray 
fluxes in order to cover as much as possible the Lx jz plane 
and hence to derive their X-ray luminosity function (LF). 

In this paper we report the results of the spectroscopic 
identifications of one of such brighter samples. The X-ray 
sources have been detected by the BeppoSAX-MEiCA in¬ 
struments in the 5-10 keV band in the framework of the 
High Energy LLarge Area Survey (HELLAS). Preliminary 
results have been presented in Fiore et al. (1999) and La 
Franca et al. (2001). The whole survey and the cata¬ 
logue is described by Fiore et al. (2001). The data have 
been analyzed in the framework of the synthesis models 
for the X-ray background by Gomastri et al. (2001), and 
the correlation with the soft X-rays has been investigated 
by Vignali et al. (2001). 

In section 2 we describe our X-ray and optical observa¬ 
tions. In section 3 we present an analysis of the evolution 
of AGN in the 2-10 keV band. Because of the reasons 


previously described, the selection and definition of type 
2/absorbed sources is still not clear, and thus we restricted 
our evolutionary studies to type 1 AGN only. The results 
are discussed in section 4. 

2. THE SPECTROSCOPIC IDENTIFICATIONS 

The spectroscopic follow up of the HELLAS sources has 
been carried out in a subsample enclosed in a region with 
6 < 79°, and outside 5^ < a < 6.5^ and 17^ < a < 20^. 
In this region the number of sources is 118 out of a total 
of 147. Their flux distribution is shown in Figure 1 and 
the sky coverage is shown in Figure 2 and listed in Table 
1 . 

The BeppoSAX X-ray positions have an uncertainty of 
about 1-1.5 arcmin, the larger at larger off-axis distances. 
We have thus searched for optical counterparts having R 
magnitude brighter than 21.0 in a circular region of 1-1.5 
arcmin of radius around the HELLAS positions (see below 
and section 3.1.1 for a discussion on the choice of this opti¬ 
cal limit). In the case of large off-axis distances, the larger 
error-boxes (1.5') have been used. 25 sources have been 
identified with cross-correlation with existing catalogues 
(NED), and 49 have been investigated at the telescope. 
The total resulting sample of 74 sources has been built up 
in such a way that to a) randomly sample the flux distri¬ 
bution of the parent catalogue of 118 sources (see Figure 
1), and b) to avoid possible biases introduced by the cross¬ 
correlation with existing catalogues (see appendix A for a 
detailed discussion). 



S (5-10 keV, erg s“‘) 


Fig. 2 .— The sky coverage of the sample of 118 sources 
from HELLAS used in this analysis. 

The follow-up observations have been carried out using 
the EFOSG2 at the 3.6m ESO telescope. Previously pub¬ 
lished identifications (Fiore et al. 1999) have been carried 
out with the RG spectrograph (RGSP) at the Kitt Peak 
4m telescope, the FAST spectrograph at the Whipple 60" 
telescope, and the Hawaii 88" telescope. Long slit and 
multi object spectroscopy have been carried with resolu¬ 
tion between 7 and lOA . The reduction process used stan- 
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Table 2 
Observing log 


Date 


Telescope 

1998 Dec 28: 

3 nights 

ESO/3.6m EFOSC2 

1999 Jun 8; 

3 nights 

ESO/3.6m EFOSC2 

2000 Jan 3; 

3 nights 

ESO/3.6m EFOSC2 

2000 Jul 29: 

3 nights 

ESO/3.6m EFOSC2 


dard MIDAS facilities (Banse et al. 1983). The raw data 
were sky-subtracted and corrected for pixel-to-pixel vari¬ 
ations by division with a suitably normalized exposure of 
the spectrum of an incandescent source (flat-field). Wave¬ 
length calibrations were carried out by comparison with 
exposures of He-Ar, He, Ar and Ne lamps. Relative flux 
calibration was carried out by observations of spectropho- 
tometric standard stars (Oke 1990). The campaign dates 
of the observations carried out at ESO are listed in Table 
2 . 

We have divided the identified sources in AGNl, AGN2, 
galaxies with narrow emission lines showing moderate-to- 
high degree of excitation (ELG), BL Lac objects, clusters 
of galaxies and stars. AGN2 includes AGNl.5, AGNl.8 
and AGNl.9. No distinction has been made between radio- 
loud and radio quite objects. The radio properties of the 
HELLAS sources will be discussed in a forthcoming paper 
(P. Giliegi, in preparation). The R-band magnitudes of 
these identifications have been pushed down to the limit 
of obtaining about less than 10% probability of having 
a spurious identification, due to chance coincidences, for 
each of the HELLAS sources (i.e. 90% reliability for the 
whole sample). Eor AGNl we have chosen a limit of R=21 
where the surface density is ~60-70/deg^ (Zitelli et al. 
1992). Tresse et al. (1996) found that 17% of galaxies 
with z < 0.3 host an emission line spectrum not typical of 
excitation from starburst (but typical of AGN2 activity), 
this percentage reduces to at minimum 8% if the maximum 
likely effect of absorption under the Balmer lines is taken 
into account. They also found that the galaxies having 
forbidden emission lines without being AGN2 are at least 
36%. These objects are mainly starburst galaxies and low 
ionization narrow emission line galaxies (LINER). We will 
call here all these sources simply emission line galaxies 
(ELG). We have assumed roughly constant these percent¬ 
ages and thus chosen a limit of R=19 for AGN2 where the 
surface density of all galaxies is about 500/deg^ and con¬ 
sequently the expected density of AGN2 is ~70-80/deg^. 
Eor ELG we have chosen a limit of R=17.5 where the sur¬ 
face density of all galaxies is about 160/deg^, and therefore 
the expected density of ELG is ~60/deg^. 

The classification of the narrow emission line galaxies 
has been carried out using standard line ratio diagnos¬ 
tics (Tresse et al. 1996). Some of the spectra showed Hq, 
in emission and H^ having a) very small emission, or b) 
absent, or c) small absorption. Such spectra show how¬ 
ever spectral features typical of strong-emission-line spec¬ 
tra such as [OH]A3727, [OHI]A4959,5007, and [SH]A6725. 
In such cases the [OHI]/H^ push the diagnostic ratio in 


the locus of AGN2, but stellar absorption at H/j could be 
significant affecting the classification. Eollowing Tresse et 
al. (1996), we have tested the classification under the as¬ 
sumption of a possible absorption of EW=2A over H^. 

We have spectroscopically identified 36 out of the 49 
HELLAS sources whose field has been investigated at the 
telescopes. For 13 sources we have found no optical coun¬ 
terparts brighter than the chosen magnitude limits (see 
discussion above). The astrometric accuracy of the new 
generation X-ray telescopes (Chandra and XMM-Newton) 
would probably allow an unambiguous identification of 
the counterparts of these HELLAS sources, but in the 
framework of this work these fields have been declared 
empty. These sources have harder X-ray spectra than 
the total sample. They have an average softness ratio of 
< >= —0.5 ± 0.3 in comparison with the average 

value —0.1 ± 0.4 of the total sample of 118 sources, where 
S and H are the 1.3-4.5 keV the 4.5-10 keV counts, re¬ 
spectively (see Fiore et al. 2001). The absence of bright 
counterparts and the average harder X-ray spectrum fa¬ 
vor the hypothesis that most of these sources preferably 
harbor absorbed AGN. If these sources correspond to opti¬ 
cal ’’normal” galaxies, such as those observed by Chandra 
and XMM-Newton (see e.g. Fiore et al. 2000; Barger et 
al. 2001; Tozzi et al. 2001; Hornschemeier et al. 2001; 
Comastri et al. 2002), we would have not been able to 
properly identify them inside our error-boxes. However, 
our sample is statistically well defined, and from the num¬ 
ber of empty fields we can directly derive an upper limit 
of ~18% (13/74) for the presence of these optical normal 
galaxies in our survey. 

The list of the observed HELLAS sources with their 
most probable spectroscopic identification is shown in Ta¬ 
ble 3. 80% of our identifications are inside an error-box of 
60 arcsecs. 2 bright galaxies (R~15-16) are at a distance 
of about 110". The 13 empty fields are listed in Table 
4. The optical spectra are shown in Figure 3; as a refer¬ 
ence, a list of the most typical emission lines for AGN are 
over-plotted with the corresponding redshift. In Table 5 
the EW, FWHM, FWZI in the observed frame of the most 
relevant emission lines are listed. 

In total, 63% (74/118) of our HELLAS subsample has 
been searched for spectroscopic identification. 61 have 
been identified: 37 AGNl, 9 AGN2, 5 ELG, 6 Glusters, 
2 BL Lac, 1 Radio Galaxy and 1 Star. 
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Table 3 

Identified Hard X-ray Sources in the HELLAS Catalogue 


BeppoSAX position 
RA (2000) DEG 
h m s o / // 

h 

Optical position 
RA (2000) DEC 
ms ° ' 

II 

^XO 

n 

Off,, 

/ 

Fx 
5-10 keV 
(1) 

R 

Spectr. ID 

z 

00 

26 

36.5 

-19 

44 

14 

00 

26 

36.1 

-19 

44 

16 

7 

13 

3.4 

18.1 

AGN2 

0.238 

00 

27 

09.9 

-19 

26 

31 

00 

27 

09.8 

-19 

26 

14 

18 

6 

1.8 

17.7 

AGNl 

0.227“ 

00 

45 

49.6 

-25 

15 

14 

00 

45 

46.3 

-25 

15 

50 

59 

24 

3.3 

17.5 

AGN2 

0.111 

01 

21 

56.8 

-58 

44 

05 

01 

21 

56.9 

-58 

44 

42 

37 

16 

2.6 

16.8 

AGN2 

0.118 

01 

34 

33.3 

-29 

58 

38 

01 

34 

33.8 

-29 

58 

16 

24 

6 

0.9 

18.0 

AGNl 

2.217 

01 

35 

30.2 

-29 

51 

22 

01 

35 

32.7 

-29 

52 

02 

53 

8 

0.9 

17.7 

AGNl 

1.344 

01 

40 

08.9 

-67 

48 

13 

01 

40 

14.7 

-67 

48 

54 

53 

8 

2.8 

12.4 

Star 

0.000 

02 

42 

01.8 

+00 

00 

46 

02 

42 

00.9 

+00 

00 

22 

29 

10 

1.5 

18.6 

AGNl 

1.112t 

03 

15 

45.0 

-55 

29 

26 

03 

15 

47.5 

-55 

29 

04 

31 

15 

2.7 

17.9 

AGNl 

0.464t 

03 

17 

32.4 

-55 

20 

12 

03 

17 

32.7 

-55 

20 

26 

14 

21 

4.1 

17.5 

AGNl 

0.406t 

03 

33 

09.6 

-36 

19 

40 

03 

33 

12.2 

-36 

19 

48 

33 

15 

4.0 

17.5 

BLac 

0.308t 

03 

34 

07.4 

-36 

04 

22 

03 

34 

07.5 

-36 

04 

04 

19 

5 

1.9 

20.1 

AGNl 

0.904 

03 

36 

51.3 

-36 

15 

57 

03 

36 

54.0 

-36 

16 

06 

34 

15 

3.7 

17.7 

AGNl 

1.537t 

04 

37 

14.5 

-47 

30 

58 

04 

37 

11.8 

-47 

31 

48 

57 

16 

2.7 

17.3 

AGNl 

0.142 

04 

38 

47.9 

-47 

29 

07 

04 

38 

47.0 

-47 

28 

02 

67 

20 

4.7 

20.5 

AGNl 

1.453 

06 

46 

39.3 

-44 

15 

35 

06 

46 

37.6 

-44 

15 

34 

19 

17 

4.3 

16.6 

AGNl 

0.153 

07 

21 

29.6 

+71 

14 

04 

07 

21 

36.2 

+71 

13 

24 

52 

8 

0.8 

17.7 

AGNl 

0.232t 

07 

41 

40.3 

+74 

14 

58 

07 

41 

50.0 

+74 

14 

48 

41 

23 

30.7 


Clust. 

0.216t 

07 

43 

09.1 

+74 

29 

20 

07 

43 

12.6 

+74 

29 

36 

22 

7 

6.0 

16.4 

AGNl 

0.312t 

08 

37 

37.2 

+25 

47 

49 

08 

37 

37.1 

+25 

47 

52 

4 

12 

2.6 

16.9 

AGNl 

0.077 

08 

38 

59.9 

+26 

08 

14 

08 

38 

59.2 

+26 

08 

14 

10 

23 

16.4 

15.3 

ELG 

0.048 

10 

32 

15.8 

+50 

51 

04 

10 

32 

16.2 

+50 

51 

20 

18 

10 

3.1 

15.9 

AGNl 

0.174 

10 

54 

19.8 

+57 

25 

09 

10 

54 

21.2 

+57 

25 

44 

38 

13 

2.6 

18.5 

AGN2 

0.205t 

11 

01 

46.4 

+72 

26 

11 

11 

01 

48.8 

+72 

25 

38 

36 

22 

7.3 

16.7 

AGNl 

1.460t 

11 

02 

37.2 

+72 

46 

38 

11 

02 

36.8 

+72 

46 

40 

3 

21 

7.9 

15.1 

AGNl 

0.089t 

11 

06 

14.0 

+72 

43 

16 

11 

06 

16.1 

+72 

44 

14 

59 

9 

1.7 

18.5 

AGNl 

0.680t 

11 

18 

11.9 

+40 

28 

33 

11 

18 

13.8 

+40 

28 

38 

23 

4 

0.8 

18.7 

AGNl 

0.387 

11 

18 

46.2 

+40 

27 

39 

11 

18 

48.7 

+40 

26 

48 

59 

5 

1.4 

18.5 

AGNl 

1.129 

12 

04 

07.6 

+28 

08 

31 

12 

04 

04.0 

+28 

07 

24 

82 

16 

5.1 


Clust. 

0.167t 

12 

17 

50.3 

+30 

07 

08 

12 

17 

52.1 

+30 

07 

02 

25 

20 

3.5 

14.0 

BLLac 

0.237t 

12 

18 

55.0 

+29 

58 

12 

12 

18 

52.5 

+29 

59 

02 

59 

13 

2.0 

18.6 

AGN2 

0.176 

12 

19 

45.7 

+47 

20 

43 

12 

19 

52.3 

+47 

20 

58 

69 

8 

1.2 

19.3 

AGNl 

0.654t 

12 

22 

06.8 

+75 

26 

17 

12 

22 

07.0 

+75 

26 

18 

2 

7 

2.5 


Clust. 

0.240t 

12 

40 

26.0 

-05 

13 

20 

12 

40 

27.8 

-05 

14 

02 

50 

12 

3.1 

18.8 

AGNl 

0.300 

12 

40 

29.6 

-05 

07 

46 

12 

40 

36.4 

-05 

07 

52 

102 

17 

1.9 

15.2 

ELG 

0.008 

13 

04 

38.2 

-10 

15 

47 

13 

04 

35.6 

-10 

15 

48 

39 

6 

1.4 

20.1 

AGNl 

2.386 

13 

05 

32.3 

-10 

32 

36 

13 

05 

33.0 

-10 

33 

20 

45 

22 

19.3 

14.9 

AGNl 

0.278t 

13 

36 

34.3 

-33 

57 

48 

13 

36 

39.0 

-33 

57 

58 

60 

22 

3.2 

10.5 

RadioG. 

0.013t 

13 

42 

59.3 

+00 

01 

38 

13 

42 

56.5 

+00 

00 

58 

59 

21 

3.2 

18.7 

AGNl 

0.804t 

13 

48 

20.8 

-30 

11 

06 

13 

48 

19.5 

-30 

11 

56 

52 

14 

2.2 

15.3 

AGN2 

0.128 

13 

48 

45.4 

-30 

29 

37 

13 

48 

44.7 

-30 

29 

46 

13 

14 

5.1 

17.1 

AGNl 

0.330 

13 

50 

09.4 

-30 

19 

55 

13 

50 

15.4 

-30 

20 

10 

80 

11 

5.1 

16.5 

ELG 

0.074 

13 

53 

54.6 

+18 

20 

33 

13 

53 

54.4 

+18 

20 

16 

18 

18 

6.8 

17.3 

AGNl 

0.217 

14 

17 

12.5 

+24 

59 

29 

14 

17 

18.8 

+24 

59 

30 

86 

13 

0.7 

19.5 

AGNl 

1.057t 

14 

18 

31.1 

+25 

11 

07 

14 

18 

31.2 

+25 

10 

50 

17 

9 

6.1 


Clust. 

0.240t 

15 

19 

39.9 

+65 

35 

46 

15 

19 

33.7 

+65 

35 

58 

41 

14 

9.4 

14.4 

AGN2 

0.044 

15 

28 

47.3 

+19 

39 

10 

15 

28 

47.7 

+19 

38 

54 

18 

5 

1.6 

20.3 

AGNl 

0.657*’ 

16 

26 

59.9 

+55 

28 

21 

16 

26 

59.0 

+55 

27 

24 

57 

11 

12.1 


Clust. 

0.130t 

16 

34 

11.8 

+59 

45 

29 

16 

34 

18.5 

+59 

45 

44 

53 

3 

0.8 

19.0 

AGN2 

0.341“ 

16 

49 

57.9 

+04 

53 

33 

16 

50 

00.0 

+04 

54 

00 

42 

20 

9.5 


Clust. 

0.154t 

16 

50 

40.1 

+04 

37 

17 

16 

50 

42.7 

+04 

36 

18 

71 

25 

12.2 

14.6 

AGN2 

0.031 

16 

52 

38.0 

+02 

22 

18 

16 

52 

37.5 

+02 

22 

06 

14 

5 

0.7 

20.7 

AGNl 

0.395 

20 

42 

47.6 

-10 

38 

31 

20 

42 

53.0 

-10 

38 

26 

80 

21 

5.3 

17.9 

AGNl 

0.363 

20 

44 

34.8 

-10 

27 

34 

20 

44 

34.8 

-10 

28 

08 

35 

15 

2.0 

17.7 

AGNl 

2.755 

22 

26 

30.3 

+21 

11 

57 

22 

26 

31.5 

+21 

11 

38 

26 

14 

3.9 

17.6 

AGNl 

0.260 

23 

15 

36.4 

-59 

03 

40 

23 

15 

46.8 

-59 

03 

14 

85 

8 

2.0 

11.2 

ELG 

0.044t 

23 

19 

22.1 

-42 

41 

50 

23 

19 

32.0 

-42 

42 

28 

116 

21 

5.7 

16.5 

ELG'* 

0.101 

23 

27 

28.7 

+08 

49 

30 

23 

27 

28.7 

+08 

49 

26 

4 

7 

0.5 

18.5 

AGN1= 

0.154 

23 

29 

02.4 

+08 

34 

39 

23 

29 

05.8 

+08 

34 

16 

56 

20 

2.9 

20.3 

AGNl 

0.953 

23 

31 

55.6 

+19 

38 

34 

23 

31 

54.3 

+19 

38 

36 

19 

17 

3.7 

18.8 

AGNl 

0.475 

23 

55 

53.3 

+28 

36 

06 

23 

55 

54.3 

+28 

35 

58 

16 

12 

4.2 

17.9 

AGNl 

0.731t 


Note.— (1) 10“^® ergcm“^ t From cross-correlation with existing catalogues; a) From Colafrancesco 
et al. (2000); b) From Gorosobel et al. (1998), misidentified in Fiore et al. 2000; c) Wrong coordinates and 
magnitudes in Fiore et al. (1999); d) Also possible AGNl: has a rest frame FWHM~1000 Km/s; e) 

__ u . ^ r _— ^ /r * i can xy — 



HELLAS -V: THE X-RAY POPULATION AND EVOLUTION 


5 


Table 4 

Empty fields in the HELLAS Catalogue 


BeppoSAX position 


RA (2000) DEC 

5-10 keV 

h m s o ! n 

(10~^® erg cm“^ 


01 

34 

28.6 

1 

CO 

O 

06 

04 

1.33 

01 

34 

49.6 

-30 

02 

34 

0.71 

09 

46 

05.3 

-14 

02 

59 

2.82 

09 

46 

32.8 

-14 

06 

16 

3.22 

13 

04 

24.3 

-10 

23 

53 

1.28 

13 

04 

45.1 

-05 

33 

37 

1.26 

13 

48 

24.3 

-30 

25 

47 

3.15 

14 

11 

58.7 

-03 

07 

02 

3.93 

22 

03 

00.5 

-32 

04 

18 

2.83 

22 

31 

49.6 

-til 

32 

08 

1.94 

23 

02 

30.1 

-1-08 

37 

06 

2.67 

23 

02 

36.2 

-1-08 

56 

42 

3.17 

23 

16 

09.8 

-59 

11 

24 

1.32 


Note.— Within 1 arcmin there are no AGNl 
down to R=21, no AGN2 down to R=19 and no 
ELG down to R=17.5. 



i°§(h5l(ev/hR) 

Fig. 4.— The histogram of the Log(L 5 keVI^ r) ratio 
for a sample of 50 (25 X-ray selected and 25 optically se¬ 
lected) AGNl (see text). We fitted the distribution with a 
Gaussian having mean -4.58 and standard deviation 0.61. 
This distribution has been assumed as representative of 
the whole population of type 1 AGNs at any redshift and 
luminosity, and used to compute the fraction of AGNl 
fainter that R=21 missed from the optical spectroscopic 
identifications of the HELLAS sample. 


3. ANALYSIS ON THE EVOLUTION OF AGNl 

Because of the difficulties previously described in quan¬ 
tifying all the selection biases in building up samples of 
type 2 (and/or absorbed) AGN, we have decided to limit 
our evolutionary analysis to AGNl only. See appendix A 
for a discussion on the absence of substantial biases on our 
sample of AGNl. We have assumed Ho=50 km/s/Mpc and 
either the (Hrra,f^A)=(l-0,0.0) or the (Um,UA)=(0.3,0.7) 
cosmologies. We have computed the 2-10 keV luminos¬ 
ity for AGNl assuming a typical spectral energy distribu¬ 
tion as computed by Pompilio, La Eranca, & Matt (2000), 
which assumes a power slope in energy with index a=0.9 
(^ (X E~°‘) and takes into account the reflection. This 
spectrum is roughly approximated by a single power a=0.6 
in the range 2-50 keV (i.e. up to z=4). This is in agree¬ 
ment with the average slope of the spectra of AGNl from 
the HELLAS sample (Vignali et al. 2001). 

3.1. The data 

3.1.1. The completeness of the sample of AGNl from 

HELLAS 

As our spectroscopic identifications are limited to mag¬ 
nitudes brighter than i? = 21, in order to estimate the den¬ 
sity of AGNl, we have evaluated which is their complete¬ 
ness in our sample. We have assumed that all AGNl follow 
the same distribution of the ratio of the 5 keV to R-band 
optical luminosity Log(L 5 keV/^ r) without dependencies 
on luminosity and/or redshift. We have estimated the re¬ 
lationship between optical and hard X-ray luminosity for 
AGNl as follow. We selected a mixed sample of both 25 
optical and 25 hard X-ray selected type 1 AGN, mostly 
in the redshift range 0<z<2, from Mineo et al. (2000), 
George et al. (2000) and Akiyama et al. (2000). The 25 
optically selected AGNl have an average Log(L 5 keV/^ r) 
ratio of -4.78 with a standard deviation of 0.59, while the 
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25 X-ray selected AGNl have an average Log(L 5 hevI^ b) 
ratio of -4.38 with a standard deviation of 0.58. In Figure 4 
the histogram of the ratio of 5 keV and R-band optical lu¬ 
minosity of the total sample of 50 AGNl at z=0 is shown. 
We fitted the distribution with a Gaussian having mean - 
4.58 and standard deviation 0.61. We have estimated that 
with this distribution, a limit of i? = 21 corresponds to 
only 1 AGNl lost from our HELLAS survey, and conse¬ 
quently we can say that a limit of i? = 21 assure a high 
level of completeness for the identification of AGNl at the 
fluxes of our HELLAS survey. At variance, as already dis¬ 
cussed, our HELLAS identification programme is biased 
against absorbed type 2 sources which probably harbor in 
the 13 empty fields found. 

Although the incompleteness for AGNl in the HELLAS 
sample is quite low, during the computation of the LF 
of AGNl, we have estimated the correct area coverage for 
our survey by multiplying the area coverage corresponding 
to each point of the Lx/z plane for: a) the spectroscopic 
completeness of the HELLAS sample (74/118), and b) the 
fraction of found AGNl as a function of z, as described 
before. 

3.1.2. The other samples 

In order to increase the statistical significance of our 
analysis, the HELLAS data have been combined to other 
hard X-ray samples of AGNl identihed by Grossan (1992), 
Boyle et al. (1998), and Akiyama et al. (2000). 

The sample of Grossan (1992) consists of 84 AGNl and 
12 AGN2 detected by HEAOl, predominantly at low red- 
shift z < 0.3. The sample covers an area of 26919 deg^ 
down to a flux limit of 1.8 x 10“^^ erg cm“^ s“^. The 
quoted fluxes at 5 keV have been converted to fluxes in 
the 2-10 keV assuming a = 0.6. 

The sample of Boyle et al. (1998) consists of 12 AGNl 
and 6 AGN2 detected by ASCA. 2-10 keV counts have 
been converted to fluxes with a count-rate-to-flux conver¬ 
sion factor of 5.8 x 10“^^ erg cm“^ s“^/count, assuming 
a = 0.6 (Georgantopoulos et ah, 1997). The sky cover¬ 
age has been taken from Georgantopoulos et al. (1997, see 
their hgure 1). 

The sample of Akiyama et al. (2000) consists of 25 
AGNl, 5 AGN2, 2 clusters, 1 star and 1 unidentified 
source. The sources have been detected by ASCA. In or¬ 
der to convert the counts of the sky coverage presented in 
their Table 1, a conversion has been applied assuming that 
1 count corresponds to 0.9 x 10“^^ erg cm“^ s“^in the 2-10 
keV band for an X-ray source with a slope q;=0.6. 

A total of 158 AGNl have been used in the redshift 
range 0<z<3. Their luminosity-redshift distribution is 
shown in Figure 5. In this distribution a gap between 
the sample of Grossan from HEAOl and the other fainter 
samples from ASCA and BeppoSAX is evident. This is be¬ 
cause surveys at fluxes in the range 10“^^-10“^^ erg cm“^ 
s“^in the 2-10 keV energy band where some hundreds of 
deg^ need to be covered are still missing. These surveys 
are in program with XMM-Newton (see e.g. Barcons et al. 
2001 ). 

3.2. The method 

The best-ht parameters for the 2-10 keV LF and its cos¬ 
mological evolution have been derived by minimization of 


the computed by comparison of the observed and ex¬ 
pected number of AGNl in 6x3 bins in the Lxjz plane. 
The adopted luminosity/redshift grid was: 6 LogLx bins, 
with ALogLx = l in the range 42.2-48.2, and 3 redshift bins 
(0.0-0.2, 0.2-1.0, 1.0-3.0). Table 6 summarizes the results 
for a number of different models. The errors quoted for 
the parameters are 68% (Ict) conhdence intervals. They 
correspond to variations of Ax^=1.0, obtained perturbing 
each parameter in turn with respect to its best-fit value, 
and looking for a minimization with the remaining param¬ 
eters free to float. The fitted LF have been tested with 
the 2D Kolmogorov-Smirnov (2DKS) test. Two different 
functional forms have been used. Following Geballos and 
Barcons (1996) and Boyle et al. (1998), we used the pure 
luminosity evolution (PLE) for the AGNl LF. We also 
used the luminosity dependent density evolution (LDDE) 
model similar to the one htted in the soft X-rays by Miyaii 
et al. (2000). 

The local (^=0) AGNl LF used for the PLE model has 
been represented by a two-power-law: 

:Lx<ii(^ = 0), 

:Lx>ii(^ = 0), 

where Lx is expressed in units of lO'^'^ergs”^. A standard 
power-law luminosity evolution model was used to param¬ 
eterize the cosmological evolution of this LF: L*^{z) = 
L*^mixz)K 



42 43 44 45 46 47 48 


LogL^(erg 2-10 keV) 

Fig. 5.— The luminosity/redshift distribution for 
all the 158 AGNl used to compute the hard X-ray 2- 
10 keV luminosity function. Ho=50 km/s/Mpc and 
(Dm,flA)=(L0,0.0) have been assumed. It includes 37 
AGNl from HELLAS, 84 AGNl from Grossan (1992), 
12 from Boyle et al. (1998) and 25 from Akiyama et al. 
( 2000 ). 


In the case of the LDDE model, as an analytical expres¬ 
sion of the present day (z=0) luminosity function, we used 
the smoothly connected two power-law form: 
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d^{L^, z = 0) 
dLogLx 


A[{L^/L,y^ + {L^/L^yT^. 


The description of the evolution is given by a factor 
e(Lx, z) such that: 

d^{L^,z) _ d^{L^,z = 0) 
dLogLx dLogLx ^ x, z , 

where 
e{L^,z) = 

( (1 -f 2)max(0,pl-aLog [LJL^]) ^ ^va) 

il + z)P^ {z<Zc-,L^>Ly 

{ e{L^,Zc)[{l +z)/{l +z^yf^ {z>Zc). 

The parameters are defined as in Miyaji et al. (2000). In 
particular, a is proportional to the decrease of the density 
evolution of the faint objects. In this parameterization of 
the LF Lx is expressed in units of ergs“^. In both mod¬ 
els the LF has been computed for luminosities brighter 
than LogL2_iofceV=42.2. Using this parameterization, for 
(Dm,U a)=( 1 . 0 , 0 . 0 ), a fit was done by keeping fixed the 
parameters Zc=1.55 (the redshift at which the evolution 
stops), a=2.5, and the ratio La/T* = 10° "^^ to the value 
found by Miyaji et al. (2000) for AGNl only (their LDDEl 
model) in the soft X-rays, and leaving all the remaining 
parameters (the two slopes 71 and 72 , the break luminosity 
L*, and the speed of the density evolution pi) free to vary. 
In the case of an (Om,DA)=(0.3,0.7) Universe, as the fit for 
AGNl only from Miyaji et al. (2000) was not available, 
their parameters from the fit for all AGNs (AGN1-I-AGN2) 
have been used: Zc=1.58, q;= 2.6, La/T* = 10° ®®. 



42 43 44 45 46 47 48 

LogL,,(erg s“‘: 2-10 keV) 


Fig. 6 .— The luminosity function of the 37 HELLAS 
AGNl. As a reference, the dotted lines show the previ¬ 
ously estimated luminosity function in the 2-10 keV from 
Boyle et al. (1998). 


3.3. Results 

3.3.1. The LF in the (flm,^A)=(L0,0.0) universe 


We started our computation of the LF in an 
(Dm,U a)=( 1.0,0.0) cosmology. Previous estimates of the 
shape and evolution of the LF of AGNl in the 2-10 keV 
range are from Ceballos and Barcons (1996), and Boyle et 
al. (1998). Boyle et al. (1998) combined the local sample 
of 84 AGNl observed by HEA Ol from Grossan (1992) with 
a fainter sample of 12 AGNl observed by ASCA (see their 
distribution in the Lx-z plane in Figure 4). They found 
that the 2-10 keV AGN X-ray LF is best represented by 
a two-power-law function evolving according to a pure lu¬ 
minosity evolution (PLE) model: Lx{z)(x{l + (see 

model 0 in Table 6 ). 

In Figure 6 the LF from only the 37 AGNl from HEL¬ 
LAS in three redshift intervals ( 0 . 0 < 2 < 0 . 2 , 0 . 2 < 2 < 1 . 0 , 
and 1.0<z<3.0) is shown. For the sake of comparison, 
the best-fit LF from Boyle et al. (1998) is also shown. 
The data have been represented by correcting for evolu¬ 
tion within the redshift bins as explained in La Franca & 
Gristiani (1997). The error bars in the figures are based 
on Poisson statistics at the 68 % confidence level. 
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42 43 44 45 46 47 48 

LogL,,(erg s^*: 2-10 keV) 

Fig. 7.— The LF of the total sample of 158 AGNl includ¬ 
ing the 37 AGNl from HELLAS, 84 from Grossan (1992), 
12 from Boyle et al. (1998) and 25 from Akiyama et al. 
(2000). The excess of faint AGNl at high redshift is even 
more evident than in the previous figure. The dotted line 
is the PLE fit by Boyle et al. (1998). The dashed line is 
our best PLE fit (see Table 6 ). 


' I ' I ' I 
z = 0.1 z==0.5 z-1.5 

' ■ V ■■ - 

A' I 
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S A 


-PLE (this work) \ 

. Boyle et al. (1998) 

H„=50, (n„,np={i.0,0.0) 

_I_,_I_,_I_,_I_,_1_ 


The data are in rough agreement with the previous es¬ 
timate of Boyle et al. (1998), but show an excess of faint 
AGNl at redshift larger than z~1.0. This feature is still 
present after combining our data with the other AGNl 
samples from Grossan (1992), Boyle et al. (1998), and 
Akyiama et al. (2000) (Figure 7). 
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Table 6 

The 2-10 keV AGN LF parameters 


Model 


71 

72 

LogL* 

k or pi 

^cut 


Prob 

fxRB^ 

0) PLE from Boyle 

1 . 0 , 0.0 

1.73 

2.96 

44.16 

2.00 


8.2 X 10“'^ 

0.04® 

24 

1) PLE 

1 . 0 , 0.0 

1.83 

3.00 

44.17 

2.12 


9.8 X 10“'^ 

0.22 

35 

2) PLE with Zcut 

1 . 0 , 0.0 

1.87 

3.03 

44.17 

2.52 

1.39 

8.9 X 10“'^ 

0.31 

34 

3) LDDE 

1 . 0 , 0.0 

0.68 

2.02 

44.03 

4.66 

1 . 55 = 

2.0 X 10-® 

0.47 

40 

4) PLE 

0 . 3 , 0.7 

1.93 

2.97 

44.24 

2.22 


8.7 X 10“'^ 

0.47 

56 

5) PLE with Zcut 

0 . 3 , 0.7 

1.94 

2.97 

44.24 

2.26 

2.38 

8.7 X 10“'^ 

0.70 

54 

6) LDDE 

0 . 3 , 0.7 

0.69 

1.95 

43.96 

4.56 

1.58= 

2.3 X 10“® 

0.77 

60 

7) LDDE softX 

1 . 0 , 0.0 

0.62= 

2.25= 

44.13 

5.40= 

1.55= 

1.4 X 10“® =’‘^ 

0.11 

48 

68% confidence errors 


+0.12 

-0.16 

+0.11 

-0.09 

+0.13 

-0.15 

+0.13 

-0.14 

+0.56 

-0.25 

8% 




Note.— (1) % of the XRB: I 2-10 = 1-95 x 10 erg cm ^ s ^ deg ^ (Chen, Fabian and Gendrau 1997) 
(a) Mpc“®(10^^erg 8“!)“^; (b) probability; (c) fixed; (d) Mpc“®(erg 
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-6 h 


-10 h 


-12 


42 


z = 0.5 


■ A 
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-LDDE 


LDDE from soft X-rays 
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La 

4 


H„=50, (n„,nd={i.0,0.0) 

J_^^_L 


43 44 45 46 47 

LogL,,(erg s“‘: 2-10 keV) 


48 


Fig. 8 .— The luminosity function of AGNl fitted with the 
LDDE model (dashed line). The soft X-ray LF (Miyaji et 
al. 2000) is also shown (dotted line), with the assumption 
LogL2_iofeeF=LogLo.5_2fceV+0-33 (see text). 


ization (the normalization has a Poisson uncertainty of 8%. 
See model 1 in Table 6). 
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42 43 44 45 46 47 48 

LogL,,(erg s“‘: 2-10 keV) 

Fig. 9.— The luminosity function in the 
(Dm,tlA) = (0.3,0.7) universe htted by our PLE model (see 
Table 6). 


z = 0.1 z = 0.5 z=1.5 



Ho= 50. (n^,nj=(o,3,o,7) -\ 


As already discussed, these samples all together collect 
158 AGNl. These data have a probability of 0.04 to be 
drawn from a PLE model such as that computed by Boyle 
et al. (1998) (model 0 in Table 6). We did not use the 
2DKS as this test is not appropriate in this case: it uses 
the cumulative distributions in the Lx / 2 : plane regardless 
of the normalization of the LF. Our best-fit to the data 
with a PLE model found a LF with similar slopes and 
break luminosity as that one of Boyle et al. (1998) but 
with a slightly larger evolution {k = 2 . 12 ^ 0 X 4 instead of 
fc= 2 . 00 ^ 0 ( 22)5 ^ signihcantly 20 (^ 5 )% larger normal- 


Our larger values of the evolution parameter k and nor¬ 
malization A are originated by the necessity of better ht- 
ting the observed higher density of faint AGNl at high red- 
shift. The 2DKS test gives a probability of 0.22 for this ht 
(see Table 6). A even better probability of 0.31 is obtained 
if a stop in the evolution is applied at redshift ZcMt=1-39 
(model 2 in Table 6) and a larger evolution (fc=2.52) is 
used. 
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Fig. 10.— The luminosity function in the 

(0„,0;)j)=(0.3,0.7) universe fitted by the LDDE model 
(see text and Table 6). 


at redshifts larger than 2.4. However our data contain not 
enough AGNl at redshift larger than 2 in order to obtain 
an accurate measure of the Zcut parameter (see Figure 4), 
therefore the errors in this parameter are quite large. Also 
the LDDE model obtain a satisfactory fit of the data (see 
Figure 10). 



^^,(erg s-‘ cm-q 


Although our hts are already statistically adequate, our 
PLE models do not fully describe an over-density of faint 
AGNl which is observed at high redshift. The data are 
not sufficiently faint to properly probe this part of the LF, 
but as this feature is similar to what observed in the soft 
X-rays (Miyaji et al. 2000), we tried to obtain a even 
better fit of the data, by using the luminosity dependent 
density evolution (LDDE) model similar to the one fitted 
in the soft X-rays, as described in the previous section. 
This model differentiates from the PLE especially in the 
faint part of the LF, at luminosities lower than L*. In this 
part of the LF the density evolution decreases in propor¬ 
tion to the faintness of the objects. As our data just start 
to probe the part of the LF which is fainter than L*, we 
cannot expect to hnd directly a fit to all the parameters of 
the LDDE model to our data, and we limited our analysis 
to a check of the compatibility (with limited changes) of 
the LDDE parameters found by Miyaji et al. (2000) to 
our data. With this model we found the best-fit to our 
data (model 3: LDDE, Figure 8) by keeping fixed the pa¬ 
rameters describing the stop in the evolution (zc=l-55), 
and the dependence on luminosity of the density evolution 
(a=2.5) to the value found by Miyaji et al. (2000) for 
AGNl only, and leaving all the remaining parameters free 
to vary (see section 3.2). The 2DKS test gives, for this 
model, a probability of 0.47. 

3.3.2. The (Vlm,^K)=(0.3,0.7) universe 

If an (Hm,UA)=(0.3,0.7) cosmology is assumed, the PLE 
models provide an ever better representation of the data 
in comparison with what found in the (Dm,UA)=(l.0,0.0) 
Universe (see Figure 9). The 2DKS probability is 0.70 and 
0.47, with and without the introduction of the Zcut param¬ 
eter, respectively. In this case even the simple PLE model 
obtain a quite good fit of the data, and the introduction of 
the Zcut parameter is necessary to stop the evolution only 


Fig. 11 .— The differential counts of the type 1 AGN used 
in the estimate of the luminosity function. The models are 
explained in Table 6. 


4. DISCUSSION AND CONCLUSIONS 

Thanks to the identihcation of 61 sources of the HEL¬ 
LAS sample we have been able to double the number of 
hard X-ray AGNl available for statistical analysis at fluxes 
in the range / 2 -iofceU ~ erg cm“^ s“^. In 

total we can use 74 AGNl at these fluxes (37 from HEL¬ 
LAS), which combined with the local sample of Grossan 
have allowed to show directly the shape of the LF of AGNl 
as function of redshift and measure its evolution. 

The PLE models provide satisfactory Hts of the 
data both in the (Hm,UA)=(l.0,0.0) and in the 
(Dm,UA) = (0.3,0.7) cosmologies. Our estimate of the LF 
in the (Hm,OA)=(l.0,0.0) has a significantly larger nor¬ 
malization in comparison to the previous measure from 
Boyle et al. (1998). 

The data start to probe in the hard X-rays the faint part 
of the LF where the excess of density of AGNl has been 
observed in the soft X-rays, justifying the implementation 
of the LDDE models. However, in both cosmologies, the 
statistic is not significant enough to distinguish between 
the PLE and LDDE models (see Table 6). In fact, in 
the prediction of the differential counts of AGNl shown in 
Figure 11, the models differentiates at fluxes fainter than 
f 2 -iokeV ~ 10“^^ erg cm“^ s“^, where the statistic is still 
poor. The new upcoming fainter surveys from Chandra 
and XMM-Newton will easily test which model is correct. 
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Fig. 12 .— The evolution of the density of type 1 AGN 
having LogLx(2-10 keV)>44.8 as a function of redshift. 
The soft X-ray data from Miyaji et al. (2000) have been 
over-plotted assuming a slope a = 0.6, which corresponds 
to a limit LogLx(0.5-2 keV)>44.5. The continuous line 
are the predictions of our LDDE model. 

In Table 6 the percentages of the contribution to the 2-10 
keV X-ray background are shown. The X-ray background 
has been computed integrating the LE up to z = 3.5 for 
Lx>W‘^^ ergjs. At variance with the results of 24% ob¬ 
tained from the evolution of the LE derived by Boyle et 
al. (1998), our models reproduces from ~35% up to ~60% 
of the XRB. We used / 2-10 = 1-95 x 10“^^ erg cm“^ s“^ 
deg“^ from Chen, Eabian and Gendrau (1997). The high¬ 
est percentages would probably imply that part of the ab¬ 
sorbed population necessary to reproduce the XRB is al¬ 
ready included in the AGNl at high redshift. However, 
more detailed analysis of this issue are beyond the scope 
of this paper. The percentages would decreas if a value f 
/2_io = 2.35 X 10”^^ erg cm“^ s“^ deg“^ from Vecchi et 
al. (1999) is assumed. 


It is interesting to notice that AGNl in the 2-10 keV 
range show an evolution up to z ~ 2 which is fairly well 
compatible with what observed in the soft X-rays. In fact 
we found a good Ht of the data in the (Hm,f2A)=(l.0,0.0) 
Universe if we assume exactly the same parameters found 
in the soft X-rays by Miyaji et al. (2000) for AGNl, by 
only looking for a fit with the break luminosity L* (model 
7 in Table 6). Miyaji et al. (2000) found LogL*=43.78. 
The value L* found by our Ht in the 2-10 keV band is 
LogL*=44.13. This luminosity difference implies an X-ray 
spectrum for AGNl with slope q;= 0.6, which is the same 
slope used in our computation of the LE. We note however 
that Miyaji et al. (2000) used a slope q;= 1.0 in their com¬ 
putations of the soft X-ray LE. Therefore, taken at a face 
value, the match between the soft and hard X-rays LFs im¬ 
plies that AGNl have a broad band concave spectrum get¬ 
ting steeper going toward lower energies. As already dis¬ 
cussed in the previous section, the (Um,UA)=(0.3,0.7) cos¬ 
mology for AGNl alone has not been analyzed by Miyaji 
et al. (2000). 

The agreement of the evolution measured in the soft and 
hard X-rays is also shown in Figure 12, where the evolution 
of the density of AGNl brighter than LogL2_iofeey=44.8 
is shown. The soft X-ray data from Miyaji et al. (2000) 
have been over-plotted assuming a slope a = 0.6, which 
corresponds to a limit LogLx(0.5-2 keV)>44.5. The con¬ 
tinuous line are the predictions of our LDDE model for 
(D„,Oa)=(1.0,0.0). 

Based on observations collected at the European South¬ 
ern Observatory, Ghile, ESO N°: 62.P-0783, 63.0- 
0117(A), 64.O-0595(A), 65.O-0541(A). This research has 
made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Labo¬ 
ratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
This research has been partially supported by ASI contract 
ARS-99-75, MURST grants Cofin-98-02-32, Cofin-99-034, 
Cofin-00-02-36 and a 1999 CNAA grant. 


APPENDIX 

THE BIASES INTRODUCED BY THE CROSS CORRELATIONS WITH CATALOGUES 

As already described, 25 sources have been identified by cross-correlation with existing catalogues. However, as the 
sample is not fully spectroscopic identified, the cross-correlation with the existing catalogues could alter the average 
characteristics (fluxes, percentages of classes of counterparts, etc.) of the subsample of identified sources. Namely, the 
subsample could not be representative of the whole sample. 

Let’s explain this with an example. In an ideal sample of 50 AGNl and 50 AGN2 a random identification of a subsample 
of 40 sources will identify about 50% AGNl and 50% AGN2 according to sampling errors. But if a cross-correlation is 
first made with a catalogue of only AGNl (let say 30) and later only 10 sources are randomly identified, the fraction of 
AGNl in the total identified subsample of 40 sources will be artificially increased (they will be at minimum 75%). Our 
subsample is risking the same sort of bias, and we wish to quantify it. 

A Kolmogorov-Smirnov (KS) two-sample test gives a 9% probability that the X-ray flux distribution of the 25 sources 
identified through the cross-correlation with catalogues belong to the same parent population of the total sample of 118 
sources. Although not significant, this low probability is due to the average slightly brighter fluxes of the cross-correlated 
subsample (< LogFs-io keV >= —12.4 ± 0.4) in comparison to the total sample (< LogFs-io keV >= —12.6 ± 0.3). We 
have thus tried to populate the faintest bins during the observing runs at the telescope, and, indeed, the KS test gives 
a 66% probability that the whole subsample of the 74 identified sources (included the empty fields) belong to the same 
parent population of the total sample of 118 sources (see Eigure 1). 

In this way, our observing runs at the telescope have recovered the possible alteration on the fraction of classes of 
sources which are identified. As the catalogues used for the cross correlation are mainly populated by AGNl (we used 
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NED), their fraction in our sample, as explained before, could be artificially increased. AGNl are 14 out of the 25 cross- 
correlated sources (56%). The fraction of AGNl is, as expected (but not significantly), lower for the sources observed at 
the telescope: 23 out of the 49 (47%). The two values are not statistically distinguishable from the observed fraction 50% 
of AGNl in the whole sample, which we thus consider representative of the whole HELLAS sample. 
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Table 1 
Sky coverage 


Log(Flux) 5-10 keV 
erg cm~^ 

Area 

deg^ 

-10.34 

55.51 

-11.85 

55.51 

-11.89 

55.48 

-11.93 

55.39 

-11.97 

55.28 

-12.00 

55.05 

-12.04 

54.66 

-12.08 

54.04 

-12.11 

53.16 

-12.15 

51.97 

-12.19 

50.63 

-12.22 

49.04 

-12.26 

47.13 

-12.30 

44.88 

-12.34 

42.44 

-12.37 

39.93 

-12.41 

37.49 

-12.45 

34.88 

-12.48 

31.93 

-12.52 

28.83 

-12.56 

25.67 

-12.59 

22.75 

-12.63 

20.03 

-12.67 

17.53 

-12.70 

15.23 

-12.74 

13.16 

-12.78 

11.32 

-12.82 

9.73 

-12.85 

8.30 

-12.89 

6.96 

-12.93 

5.67 

-12.96 

4.56 

-13.00 

3.60 

-13.04 

2.88 

-13.07 

2.35 

-13.11 

1.96 

-13.15 

1.60 

-13.19 

1.27 

-13.22 

0.94 

-13.26 

0.68 

-13.30 

0.47 

-13.33 

0.33 

-13.37 

0.22 

-13.41 

0.14 

-13.44 

0.08 

-13.48 

0.05 

-13.52 

0.03 

-13.56 

0.01 

-13.59 

0.00 
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Table 5 

Emission Line Measurements 


Name 

z 

Type 

Mgll 

(A) 

OH 

(A) 

H,3 

(A) 

OHI 

(A) 

H„+NH 

(A) 

SII 

(A) 

H002636-194416 

0.238 

FWZI 


62.1 

38.1 

62.1 





FWHM 


19.7 

15.4 

18.5 





EW 


82.4 

14.6 

133.1 



H004546-251550 

0.111 

FWZI 




40.1 

104.2 

44.1 



FWHM 




14.6 

69.7 

50.4 



EW 


-1.3 

-0.2 

5.6 

22.2 

2.5 

H012157-584442 

0.118 

FWZI 


44.1 

38.1 

64.1 

120.2 

48.1 



FWHM 


19.7 

13.1 

16.8 

45.0 

80.0 



EW 


30.4 

2.9 

21.4 

53.5 

20.8 

H013434-295816 

2.217 

FWZI 

“452.7 








FWHM 

“123.6 








EW 

“156.5 






H013533-295202 

1.344 

FWZI 

282.5 








FWHM 

67.6 








EW 

61.2 






H033408-360403 

0.904 

FWZI 

412.6 

40.1 







FWHM 

174.0 

26.5 







EW 

120.5 

6.1 





H043712-473148 

0.142 

FWZI 


36.0 

122.2 

64.1 

162.3 

70.1 



FWHM 


16.4 

36.4 

20.9 

39.7 

73.0 



EW 


5.2 

65.4 

23.6 

241.2 

17.0 

H043847-472802 

1.453 

372.6 









64.9 









61.4 







H064638-441534 

0.153 

FWZI 


60.1 

128.2 

60.1 

122.2 

46.1 



FWHM 


18.3 

47.2 

21.3 

36.1 

29.0 



EW 


7.9 

61.0 

78.9 

168.4 

10.0 

H083737-I-254752 

0.077 

FWZI 


52.1 

69.4 

56.8 

327.2 

60.1 



FWHM 


40.3 

28.8 

15.2 

71.8 

43.0 



EW 


51.8 

16.0 

11.3 

291.6 

17.0 

H083859-I-260814 

0.048 

FWZI 



36.7 

40.1 

110.2 

70.1 



FWHM 



11.7 

23.2 

42.1 

32.2 



EW 


-8.9 

1.9 

4.1 

46.0 

6.7 

H103216-I-505120 

0.174 

FWZI 


30.0 

203.7 

36.7 

150.2 

33.4 



FWHM 


12.0 

68.1 

13.2 

79.9 

15.0 



EW 


5.8 

65.1 

13.1 

163.2 

6.5 

H111814-I-402838 

0.387 

FWZI 

340.5 

50.1 


40.1 





FWHM 

184.0 

46.7 


10.3 





EW 

393.8 

6.7 

-0.7 

6.7 



H111849-I-402648 

1.129 

FWZI 

“200.3 








FWHM 

“65.6 








EW 

“59.9 






H121853-I-295902 

0.176 

FWZI 


43.4 

36.1 

56.8 

153.6 

50.1 



FWHM 


12.4 

12.1 

12.5 

26.8 

20.9 



EW 


25.4 

17.8 

199.6 

186.1 

22.6 

H124028-051402 

0.300 

FWZI 


80.1 

166.9 

63.4 





FWHM 


25.6 

61.3 

23.9 





EW 


36.1 

90.8 

67.1 




Note.— Observed frame. Negative values correspond to absorption lines; (a) C IV line. 
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Table 5 

Emission Line Measurements 


Name 

Z 

Type 

Mgll 

(A) 

OH 

(A) 

H/3 

(A) 

OHI 

(A) 

Hc-l-NII 

(A) 

SH 

(A) 

HI24036-050752 

0.008 

FWZI 



20.0 

-10.0 

90.1 

76.8 



FWHM 



13.9 

-7.0 

34.0 

26.0 



EW 



0.8 

-0.4 

26.3 

6.8 

HI30436-101549 

2.386 

FWZI 

“468.7 








FWHM 

“128.4 








EW 

“537.0 






HI34820-301156 

0.128 

FWZI 


52.1 


56.1 

112.2 




FWHM 


21.5 


22.0 

53.9 




EW 


19.0 

-0.4 

12.7 

63.0 


HI34845-302946 

0.330 

FWZI 


36.0 

220.4 

72.1 





FWHM 


15.8 

83.0 

20.5 





EW 


2.3 

84.2 

55.8 



HI35015-302010 

0.074 

FWZI 


60.1 


52.1 

168.3 

56.1 



FWHM 


23.1 


24.5 

64.4 

31.7 



EW 


100.5 

-0.6 

16.0 

48.8 

8.0 

HI35354+182016 

0.217 

FWZI 


52.1 

24.0 

64.1 

252.4 

40.1 



FWHM 


13.7 

9.3 

19.7 

95.0 

47.9 



EW 


47.3 

2.1 

22.8 

154.0 

12.4 

HI51934+653558 

0.044 

FWZI 


36.0 


44.1 

92.2 

44.1 



FWHM 


13.6 


9.7 

9.0 

14.5 



EW 


43.6 

-1.3 

22.6 

26.4 

6.1 

HI63419+594504 

0.341 

FWZI 

144.2 

48.1 

72.1 

56.1 





FWHM 

19.9 

12.9 

12.9 

13.7 





EW 

33.1 

19.7 

30.8 

267.2 



HI65043+043618 

0.031 

FWZI 


48.1 

36.1 

52.1 

124.2 

56.1 



FWHM 


40.1 

27.2 

21.0 

49.1 

32.6 



EW 


45.9 

2.8 

26.6 

43.0 

12.3 

HI65238+022206 

0.395 

FWZI 

184.3 

56.1 

164.3 

64.1 





FWHM 

70.0 

24.1 

129.5 

21.5 





EW 

66.6 

22.5 

37.5 

68.5 



H204253-103826 

0.363 

FWZI 


72.1 

184.3 

76.1 





FWHM 


22.5 

48.5 

22.3 





EW 


15.4 

73.0 

201.8 



H204435-102808 

2.755 

FWZI 

“248.4 








FWHM 

“101.9 








EW 

“56.6 






H222632+211138 

0.260 

FWZI 



125.5 

26.7 





FWHM 



41.6 

7.9 





EW 



36.4 

9.3 



H231932-424228 

0.101 

FWZI 


68.1 

40.1 

56.1 

172.3 

64.1 



FWHM 


23.1 

26.3 

21.8 

46.9 

38.4 



EW 


20.4 

5.0 

24.2 

77.6 

7.5 

H232729+084926 

0.154 

FWZI 


28.0 

92.2 

60.1 

116.2 

56.1 



FWHM 


15.9 

80.3 

18.7 

34.2 

37.5 



EW 


13.0 

16.8 

19.1 

60.9 

16.2 

H232906+083416 

0.953 

FWZI 

368.6 








FWHM 

76.1 








EW 

133.6 






H233154+193836 

0.475 

FWZI 

344.5 








FWHM 

75.0 








EW 

93.5 
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Fig. 3.— The optical spectra of the most probable identified counterparts of the HELLAS sources, 
emission lines for AGN are over-plotted with the corresponding redshift. 
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Fig. 3.— continued 








































































































